Electronic structure of graphene on single crystal copper substrates 
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The electronic structure of graphene on Cu(lll) and Cu(100) single crystals is investigated using low energy 
electron microscopy, low energy electron diffraction and angle resolved photoemission spectroscopy. On both 
substrates the graphene is rotationally disordered and interactions between the graphene and substrate lead to 
a shift in the Dirac crossing of ~ -0.3 eV and the opening of a ^ 250 meV gap. Exposure of the samples 
to air resulted in intercalation of oxygen under the graphene on Cu(100), which formed a (v 7 ^ x 2v / 2)R45° 
superstructure. The effect of this intercalation on the graphene ty bands is to increase the offset of the Dirac 
crossing (~ -0.6 eV) and enlarge the gap (~ 350 meV). No such effect is observed for the graphene on Cu(l 1 1) 
sample, with the surface state at F not showing the gap associated with a surface superstructure. The graphene 
film is found to protect the surface state from air exposure, with no change in the effective mass observed, as for 
1 monolayer of Ag on Cu(l 1 1). 



I. INTRODUCTION 

Graphene 's unique electronic properties 1 have excited in- 
tense research interest since its isolation by micromechani- 
cal cleavage in 2004 2 . However this technique, used in these 
early investigations, is unsuitable for large-scale production 
of graphene sheets. Thermal decompositon of SiCpE, sur- 
face segregation of C dissolved in metal substra te and de- 
composition of hydrocarbons at metal surfacePM have all 
been shown to produce graphene and are potentially scalable 
to commercial production. 

Recently the growth of large scale graphene films (~ 1 
m wide) has been accomplished using chemical vapor de- 
position (CVD) onto poly crystalline Cu foils 12 , leading to 
a strong interest in the electronic and structural properties 
of the graphene/Cu interface. The structural properties of 
graphene/Cu was investigated 13 in detail, however electronic 
structure data, particularly Angle Resolved Photoemission 
Spectroscopy (ARPES), on such films are hampered by the 
many randomly oriented domains in both the graphene and the 
substrate. This leads to a large collection of electronic bands 
from both the substrate and the graphene, greatly confusing 
the interpretation of the ARPES data. 

A recent attempt to overcome this problem by Varykhalov 
et al. involved measuring graphene on 1 monolayer of 
Cu/Ni(lll), which indicated that the interaction between 
graphene and the substrate was significantly underestimated 
by current DFT L(S)DA theory 14 . The method involved grow- 
ing a graphene film on a 15-20 monolayer Ni(lll) film on 
W(110) followed by intercalating a monolayer of Cu(lll) 
underneath the graphene. While providing insight into the 
electronic structure of graphene on Cu, it is an open ques- 
tion whether this arrangement reflects the properties of the 
interface between graphene and bulk copper. The growth of 
graphene on single crystal Cu, undertaken in the current study, 
provides an ideal system for investigating the electronic struc- 



ture of graphene on Cu by eliminating the rotational disorder 
in the substrate. As the bulk and surface Cu bands are well 
known, in our samples it is easy to identify the graphene n 
bands. 

Graphene films were grown on Cu(lll) and Cu(100) sin- 
gle crystals and characterized using low energy electron 
microscopy (LEEM) and low energy electron diffraction 
(LEED). The electronic structure of the graphene samples was 
investigated using ARPES measurements. Here we show that, 
similar to graphene on 1ML Cu/Ni(lll) 14 , graphene on Cu 
(100) and (111) single crystals is n-doped with the Dirac 
crossing energy located -0.3 eV below the Fermi level 
E F . We find a gap of ~ 250 meV at the Dirac crossing 
of the graphene tt bands, significantly larger than the gap in 
graphene on 1ML Cu/Ni(lll) (~ 180 meV) and the LDA es- 
timate (~ 1 1 meV^m 

We also show that air exposure leads to oxygen interca- 
lation under the graphene on Cu(100), forming a (y/2 x 
2a/2)R45° superstructureand giving a larger doping (~ -0.6 
eV Dirac crossing shift) and an increased gap (~ 350 meV 
in the graphene tt bands at E^). The Cu(lll) surface state 
at T (Brillouin zone centre) does not show the gap expected 
for such an intercalted surface structure, confirming that expo- 
sure to air does not produce the same effect for graphene on 
Cu(lll). In contrast the graphene protects the surface state 
at T from air exposure, with no change in the effective mass 
observed. 



II. METHODS 

Cu(100) and Cu(lll) single crystals were initially cleaned 
by annealing at 900° C for about 12 hours in a flowing mixture 
of hydrogen and argon at atmospheric pressure followed by 
sputtering with argon ions. Subsequent annealing at 750° C 
in the LEEM produced surfaces with atomically flat terraces 
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separated by monatomic Cu steps and step bunches. Large 
bunches of steps formed on Cu(lll) at pinning sites spaced 
about 2-5 microns apart. Carbon was deposited onto the 
substrates at 850 — 900° C from a graphite rod heated by an 
electron beam. The growth was imaged in real time using 
LEEM. The deposition was stopped near the completion of the 
first graphene layer using the ability of LEEM to distin guish 
bare Cu, monolayer graphene and bilayer graphene^^. 

Selected-area LEED patterns were acquired after growth 
from regions 2 or 20 fim in diameter. The graphene-covered 
Cu crystals were removed from the ultrahigh vacuum of the 
LEEM system and exposed to air for about 3 hours before 
again being placed in vacuum and degassed at about 350° C for 
ARPES measurements. After the ARPES measurements, the 
samples were re-introduced into the LEEM for further anal- 
ysis. Both the graphene/Cu(100) and graphene on Cu(lll) 
specimens experienced an additional 3 hours exposure to air 
while the graphene/Cu(lll) was also kept in a desiccator for 
several weeks prior to analysis. The samples were briefly de- 
gassed at about 250° C before the subsequent LEEM/LEED 
measurements. The LEEM/ LEED measurements where per- 
formed at Sandia National Laboratories in Livermore, Cali- 
fornia. 

ARPES spectra were obtained at the Electronic Structure 
Factory endstation (SES-R4000 analyzer) at beamline 7 of the 
Advanced Light Source, Lawrence Berkeley National Labo- 
ratory. A photon energy of 95 eV was used giving overall 

resolutions of ~ 25 meV and ~ 0.01 A . The spatial area 
sampled by the ARPES spectra was typically 50 to 100 /mi. 
During the measurements the samples were cooled to ~ 20 K 
using a liquid He-cooled cryostat and the pressure was < 2 
x 10- 10 Torn 

The ARPES spectrum of graphene on Cu single crys- 
tals was modelled as follows. The bare tt band dispersion 
£?bare(k) was computed using a first nearest neighbour tight 
binding (TB) model based on the approach given by Saito et 
al.[00. The secular equation \H — J Ebare±(k).5'| = is solved 
for the eigenvalues £?bare±(k), with 



H 



£2p - Egap/% 7o/(*0 



7o/(*0* e 2p 



o/2 



(1) 



S = 



sof(k)* 



sof(k) 



and 



f(k) = ex.p(ik.Ri + ik.R2 + ik.Rs) 
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where 70 is the hopping potential, so is the overlap potential 
and it^it^, ^3 are the vectors denoting the position of the 3 
first nearest neighbor carbon atoms. 

A lattice constant, a = 2.46, is employed as are the fit- 
ted parameters, 7o=-3.24 eV and so=0.0425 eV, determined 
by Bostwick et al. for graphene on (6>/3 x 6<\/3)R30 o C- 



SiC(0001 J21 The fitted parameter, e 2p , is the offset of the 
Dirac energy, E& , from the Fermi level due to doping of the 
graphene by the substrate and is determined by comparison to 
the graphene on Cu ARPES measurements. The E gap param- 
eter is introduced by us to account for the introduction of an 
energy gap at E& due to the breaking of graphene 's A-B lat- 
tice symmetry. Broadening of the bare band is introduced to 
the model through the self-energy via the single-particle spec- 
tral function: 



A(E,h) = 



|ImS(E,k)| 
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The self-energy, E(22, k) is determined using the semi- 
empirical method of Bostwick et alPS where the linewidth 
of the ARPES data is used to determine the imaginary compo- 
nent of the self-energy, which is Hilbert transformed to get the 
real component. This experimental self-energy is then used 
to recreate the experimental data and a self-consistent fitting 
used to further refine the self-energy. An additional gaussian 
broadening term is added to the current model to account for 
experimental broadening in energy, AE = 25 meV, and mo- 
mentum, Ak = 10 mA _1 . Fig. [IJshows experimental data for 
graphene on (6>/3 x 6v / 3)R30° C-SiC(OOOl) (data adapted 
from Bostwick et al. 20 ) compared to a model spectral function 
fitted according to the above procedure. The fitted values, ex- 
cepting €2 P and Eg a p, used for the graphene on Cu data model 
are those obtained from this fit as the rotational disorder pro- 
hibits fitting in this case. 

For comparison to the graphene on single crystal Cu exper- 
imental data many rotationally disordered domains need to be 
modeled. This is done by producing a 3D (2 momentum and 
1 energy) spectral function for a single graphene domain, as 
shown in fig. [T] This single domain is then rotated around the 
T point (in this case 20 single domains, with evenly spaced ro- 
tation angles between ± 2.5 °) to produce the spectral function 
from 20 rotationally disordered graphene domains. The rota- 
tional disorder model is then produced by summing the spec- 
tral function from each of these rotational domains together. 
Inclusion of additional domains, or domains with larger rota- 
tions, did not appreciably change the resulting spectral func- 
tion, hence these values where employed. 



III. RESULTS AND DISCUSSION 

A. LEEM of as-grown films 

Figs. [2ja) and[3ja) show representative LEEM micrographs 
of graphene on Cu(l 1 1) and Cu(100), respectively. These im- 
ages show that graphene grows in grains separated by bound- 
aries, whose intensity can be bright or dark depending on the 
electron energy used for imaging. The bright dots in Fig. [2] 
(a) are the graphene nucleation sites. Similar growth struc- 
tures have been observed in STM measurements of graphene 
on Cu(100;P. 

Figs. |2jb) and|3jb) show representative LEED patterns in- 
tegrated over the 20 fim regions of Figs. |2ja) and [5Ja) for 
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FIG. 1. Experimental fermi surface and tt band structure obtained 
from graphene on (6^3 x 6 V3)R30° C-SiC (a i, ii and iii) for com- 
parison to a semi empirical model (b i, ii and iii). The model em- 
ploys a first nearest neighbour tight binding model for the bare (un- 
broadened) bands with the parameters (70 =-3. 24 eV, so =0.0425 eV 
and 62 P =-0.47 eV) determined from a fit to the experimental data. 
Broadening is performed using the self-energy, obtained by a fit to 
the experimental data in the V - K direction, and experimental broad- 
ening of 25 meV and 0.01 A -1 . The simulation and experiment agree 
qualitatively. 




FIG. 2. LEEM/LEED characterization of a nearly complete single- 
layer graphene film on Cu(l 1 1). (a) LEEM image with 20 /xm field of 
view. Bright dots are nucleation sites and bright lines are boundaries 
between graphene islands, (b) LEED from a 20-/im diameter area 
of (a). First-order diffraction spot of Cu marked by red arrow and 
portion of graphene arc marked by blue arrow. 



Cu(lll) and (100), respectively. The LEED pattern obtained 
from graphene on Cu(lll) consists not of six sharp spots ex- 
pected for a perfect graphene sheet, but instead is smeared 
into a nearly compete ring. This shows that graphene grows 
on Cu(lll) with a substantial amount of azimuthal disorder. 
The graphene diffraction is most intense near the first-order 
Cu(lll) spots, indicating a preference for graphene to align 
with the Cu lattice. Subsequent growth experiments estab- 
lished that the graphene not aligned with the Cu(l 1 1) substrate 
originated from islands that nucleated at the large bunches 
of Cu steps 2 ^. Reducing the bunch density leads to higher 
fractions of aligned graphene. The graphene LEED pattern is 
also markedly broadened in the radial direction, a result of the 
Cu hillocks that form when uncovered Cu sublimes during 



graphene growth 13 . The broadening of the graphene LEED 
pattern in the radial direction results from the fact that the 
graphene conforms to this roughened topography. 

The LEED pattern obtained from graphene on the Cu(100) 
crystal is similar to that from Cu(lll) crystals and (100) 
grains of copper foils. The LEED pattern consists of two sets 
of broad, six-fold arc^H. (One set of arcs is stronger in the 
example in Fig. [3jb)). The two sets of arcs arise because 
there are two symmetry-equivalent ways to put the six-fold 
graphene on the four-fold Cu(100) surface. The arcs occur be- 
cause the graphene has a large range of in-plane orientations. 
In addition, the arcs are broadened in the radial direction, a 
result of the roughened surface topology, as for Cu(lll) and 
STM measurements 21 . Reducing the integration area for the 
LEED patterns down to 2 microns does not completely elimi- 
nate the azimuthal disorder observed in the LEED patterns for 
both Cu(100) and Cu(lll). Therefore the size of a domain 
having a single in-plane orientation is smaller than the typical 
grains seen in Figs. [2] (a) and [3] (a). 



B. Angle Resolved Photoemission Spectroscopy 

Angle Resolved Photoemission Spectroscopy (ARPES) 
measurements of graphene on Cu(lll) and Cu(100) surfaces 
are presented in Figs. |4]and[6j respectively. The Fermi sur- 
face of the graphene on Cu(lll) is presented in Fig. |?Ja). 
The large nearly continuous arc (indicated by the large dashed 
line) corresponds to the Fermi surface of the rotationally dis- 
ordered graphene tt bands. Similar to the LEED pattern pre- 
sented above, an increase in the intensity in the Cu(lll) T-K 
direction is observed. This is the direction along which the 
two lattices are aligned.The other features in Fig. |4ja) cor- 
respond to cuts through the Cu(lll) Fermi surface and the 
central circular arc of the Cu(l 1 1) surface state 23 24 . Since the 
clean Cu(lll) surface state does not survive exposure to air 
and since the LEEM images show that the substrate is cov- 
ered, the surface state observed here must be localized to the 
interface between graphene and Cu(l 11). 

The experimental spectral function obtained along the 
white line in Fig. |4ja) is presented in Fig. [4jb) with the TB 
fit obtained using €2 P = -0.3 eV shown in blue (dashed line). 
The good agreement between the TB fit and the experimen- 
tal data indicates that the graphene is electron (n-type) doped 
by the Cu layer. The effect of this doping is a shift in the 
Dirac crossing to -0.3 eV, in good agreement with data ob- 
tained from graphene on 1ML Cu/Ni(l 11^-0.31 e V). 

Varykhalov et al. 14 also found that a gap was opened at the 
Dirac point ( ~ 180 meV) in the graphene on 1ML Cu/Ni(l 11) 
system. At first glance the spectral function in Fig. |4](b) also 
suggests a gap for graphene on single crystal Cu however the 
Energy Distribution Curves (EDC) and Momentum Distribu- 
tion Curves (MDC) shown in Fig. [4] (d) and (e) indicate in- 
tensity is present throughout the Dirac crossing region which 
suggest the presence of states. Analysis shows that this in- 
tensity is largely caused by the smearing of the bands due to 
rotational disorder, which provides an additional level of com- 
plexity. An accurate assessment of the size of the gap requires 
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FIG. 3. LEEM/LEED characterization of single-layer graphene film on Cu(100). (a) LEEM image with 20 /im field of view. Dark lines are 
boundaries between graphene islands, (b) LEED from 20 /am diameter area of as-grown graphene. One graphene arc marked by blue arrow, 
(c) LEED from a 2 /xm diameter area of air-exposed sample. Diffraction spots of Cu, graphene and intercalated O marked by red, blue and 
green arrows, respectively. 
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FIG. 4. Experimental Fermi surface (a) obtained from graphene on single crystal Cu(l 11). Also presented is the experimental spectral function 
(b), semi emperical model spectral function (c), energy distribution curves (d) and momentum distribution curves (e) obtained along the white 
line in (a) (Cu Brillouin zone T-K direction). The graphene ty bands calculated from a first nearest neighbour tight binding model (7o=-3.24 
eV, so=0.0425 eV and €2 P =-03 eV) are overlayed as the large dashed blue lines. A parabolic fit to the Cu(lll) surface states is also shown 
(small dashed yellow lines). All other unmarked bands are due to the bulk Cu. The graphene Fermi surface is observed as an uninterrupted 
ring due to rotationally disordered graphene domains. The preferential alignment (higher intensity) of the domains along the V - K direction is 
also observed in LEED. 



taking the rotational disorder into account, as discussed next. 

The predicted EDC obtained from the rotational disorder 
model described above, labelled Gm , with different simu- 
lated gaps are shown (blue solid lines) in Fig. [5] (a). The most 
significant change with gap size is reducing the depth of the 
minima at the Dirac crossing ( ~ -0.3 eV). The experimen- 
tal EDC curve is overlayed (red, dashed line) on three EDC 
curves (E g = 270, 250 and 230 meV). The best fit between 



experiment ( fig. [4] (b)) and theory ( fig. [4] (c)) is obtained 
when a gap of ~ 250 meV is used in the model. This is an 
appreciably larger gap than was found for graphene on 1 ML 
Cu/Ni(lll)^, Eg = 180 meV, which indicates a stronger in- 
teraction between the graphene and single crystal Cu. We find 
that the doping level and bandgap does not vary around the tt 
band arc. 

In contrast to the above analysis, a mere visual inspection 
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FIG. 5. The energy distribution curves (red dashed lines) obtained 
at the K point from the semi-empirical model for a range of sim- 
ulated band gaps (blue solid lines) are presented for (a) Cu(lll), 
together with a single band model Gm whose gap E g (Gm) is 
varied, and (b, c) for Cu(100), which is compared to a two-band 
model Gibo, Gf 00 . In (b), the gap E g (Gf 00 ) is varied for fixed 
E g (Gioo)= 250 meV, while in (c), the gap i^(Gibo) is varied for 
fixed E g (Gf 00 )=350 meV. In all panels, the light blue lines represent 
best-fit models to the data. 



of the data in Fig. 4b and 5a would suggest a much larger gap, 
around 350 to 400 meV. This demonstrates that the rotational 
smearing of the spectral function due to azimuthal disorder is 
to artificially enhance the gap size by up to 60% . 

The graphene on Cu(100) Fermi surface is shown in Fig. [6] 
(a), where the graphene ring is twice as broad as for graphene 
on Cu(lll), Fig. [4] (a). The spectral function along the white 
line in Fig. [6] (a), shown in Fig. [6] (b) indicates a shoulder 
(green, dash dot dot line) offset from the main band (blue, 
dash dot line). The TB fit (green, dash dot dot line) describes 
this shoulder well in [6] (a) and (b). Although the main peak 
in Fig. [6] can be described by the same TB fit, a higher Dirac 
crossing offset (e2p = -0.6 eV) is obtained by the fitting. The 
extra width of the Fermi surface can therefore be ascribed 
to two distinct doping levels for graphene on single crystal 
Cu(100). The origin of the two doping levels will be dis- 
cussed in the next section, however it is important to note that 
one doping level is the same as the graphene on single crys- 
tal Cu(lll) case (Gm in Fig. [4]and G^ 00 in Fig. |6J and an 
extra doping level most likely due to some difference in the 
sample(Gfoo m Fig- ®- 

Similar to graphene on single crystal Cu(lll), Fig. |4](d) 
and (e), the EDC and MDC plots in Fig. [6](d) and (e) indicate 



intensity in the region of the Dirac crossings. To determine 
the size of the gap in both doping levels ( Gfo and Gf 00 ) it 
is once again necessary to compare the EDC curves at the K 
point from the rotational disorder model and the experimental 
data. In this case the model is extended by creating the 3D 
spectral functions, as used in the graphene on Cu(lll) case, 
for each of the two doping levels indicated by the TB bare 
bands in Fig. [6] and summing them together. The EDC curves 
at the K point (blue, solid lines)are shown for different sim- 
ulated gaps in the G% (Fig JH (b)) and G&o (Fig- @(c)) 
bands, respectively. In Fig. [5 (b) the G^ o band gap is set 
to 250 meV and in Fig. [5](c) the Gf o band gap is set to 350 
meV. The overlying experimental EDC (red dashed lines) give 
the best fit for a gap in the Gf 00 band (Fig. [5](b)) of 350 meV 
and a gap in the Gf 00 band (Fig. [5](c)) of 250 meV. 

The model spectral function obtained based on the two TB 
bare bands from Fig. [6] and the band gaps E g (Gum) = 250 
meV and E g (Gi 00 ) = 350 meV is shown in Fig. pUc), and 
matches well with the experimental spectral function in Fig. 
[6] (b). A structural origin for the two doping levels, where 
two facets on the surface produce a momentum offset of two 
similarly doped graphene domains, was also considered. The 
model spectral function for such a system did not reproduce 
the shape of the experimental EDC plots in Fig. |5](b) and (c), 
and therefore is discounted. Such facetted features where also 
not observed in STM measurements of Rasool et alP2 

Summarizing the ARPES data, graphene bands Cu(lll) 
and Cu(100) (Gm and G^ o» resp.) were observed with iden- 
tical gap (£^=250 meV) and doping level (e2 P = -0.3 eV). A 
second graphene band (G^ 00 ) was observed on Cu(100), with 
larger gap (~ 350 meV) and doping level (~ 0.6 eV). This 
implies that the interaction of the graphene with the substrate 
layer is substantially stronger in the G^ o doping regions. 

The interactions between graphene and metal substrates has 
been modeled extensively using DFT L(S)DA theory^ and 
by considering van der Waals interactions using variations on 
vdW-DF theor\ 26 . Khomyakov et al. 25 , in addition to cal- 
culating the band structure, investigated the effect of varying 
the bond distance between the graphene layer and a metallic 
substrate. For graphene on Cu(lll) they found that varying 
the bond distance from 3 A to 4 A can shift the Dirac cross- 
ing from ~ -0.5 eV to ~ +0.25 eV. The close agreement of 
the Dirac crossing for graphene on single crystal Cu(lll), 
graphene on single crystal Cu(100) and graphene on 1ML 
Cu/Ni(l 11 jSI therefore indicates that the graphene-Cu bond 
distance is similar in all cases. 

On both Cu(lll) and (100) the spectral function does not 
vary appreciably around the arc of the graphene tt band. That 
is, doping level and bandgap do not change with graphenes 
in-plain orientation. Furthermore, the Dirac crossing values 
(dopings) and bandgaps are similar on the two Cu surfaces. 
So graphenes electronic structure is not very sensitive to the 
physical and electronic structures of the Cu surfaces or the 
precise alignment of the film/substrate lattices. These obser- 
vations suggest relatively weak interactions between the film 
and the Cu. But there are sufficient interactions to dope the 
graphene and open a bandgap by breaking the symmetry of 
the graphene lattice. The larger band gap on the Cu single 
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FIG. 6. Experimental Fermi surface (a) obtained from graphene on single crystal Cu(100). Also presented is the experimental spectral function 
(b), semi emperical model spectral function (c), energy distribution curves (d) and momentum distribution curves (e) obtained along the white 
line in (a) (Cu Brillouin zone T-K direction). The graphene tv bands calculated from a first nearest neighbour tight binding model (7o=-3.24 
eV, so =0.0425 eV and €2 P =-03 eV) are overlay ed as the dash dot blue and dash dot dot green lines, respectively. All other unmarked bands 
are due to the bulk Cu. The graphene Fermi surface is observed as an un-interuppted ring due to rotationally disordered graphene domains. 



crystals compared to graphene on 1ML Cu/Ni(lll) 14 may re- 
sult from greater substrate-induced symmetry breaking. This 
is in contrast to the DFT calculations, which predict only a 
small band gap (11 meV^l 



C. Air exposed films 

To gain insight into the origin of the two doping levels ob- 
served for graphene on Cu(100), the sample was re-examined 
by LEEM/LEED after the ARPES experiments. After de- 
gassing at about 250 °C, the LEEM images were indistin- 
guishable from the as-grown film. However, Fig. [3](c) reveals 
new LEED features - four weak, radial lines rotated 45° rel- 
ative to the directions of the first-order Cu spots. At the ends 
of the lines are weak spots (see green arrow) separated by y/2 
times the length of the Cu reciprocal lattice. This pattern can 
be interpreted as induced by the interaction of air with the 
graphene/Cu(100) system. Since graphene on other substrates 
survives air exposure and is cleaned easily with annealing at 
300° C, it is likely that the effect represents intercalation of an 
atmospheric species between the graphene and the copper. 

As shown by LEED and scanning tunneling microscopy 
(STM)pEO annealing oxygen adsorbed on Cu(100) at about 
300°C forms two domains with (y/2 x 2>/2)R45 symme- 
try at a saturation coverage of one half monolayer. Since our 
samples were annealed to a similar temperature after transfer- 



ring them through the air to the ARPES chamber, it is likely 
that we are observing the intercalation of oxygen between the 
graphene and the Cu(100) surface. In our case, the inter- 
calated oxygen is not well-ordered, giving lines rather than 
strong LEED spots. 

Auger electron spectroscopy (AES), not shown here, con- 
firmed that the air-exposed sample contained oxygen. Survey- 
ing the surface by examining LEED from areas 2 fim in di- 
ameter established that the intercalation was uniform on this 
length scale. However, we cannot exclude spatial variations 
in oxygen content at smaller length scales. Thus, we propose 
that the two doping levels of graphene on Cu(100) observed in 
ARPES arise from intercalated and non-intercalated regions, 
which would lead to two different levels of charge transfer to 
the graphene. The effect of the oxygen on the graphene then 
also acts to increase the size of the gap from ~ 250 meV to 
^350 meV. It is interesting that in their study of graphene on 
1ML X/Ni(lll), with X being Cu, Ag or Au, Varykhalov et 
slI™ also found that the band gap became larger with increas- 
ing doping. We see a similar increase in the band gap for the 
oxygen intercalated region, which has a hgiher doping level. 

Even lengthy air exposure did not change the LEEM images 
and LEED from the graphene /Cu(lll) specimen. However, 
AES revealed the presence of oxygen. Previous studies have 
established that adsorbed oxygen on Cu(lll) does not form 
ordered structures and changes the work function less than 15 
meV. 29 While we cannot totally exclude some oxygen interca- 
lation upon air exposure for graphene/Cu(lll), the net effect 
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-0.2 0.0 p.Z 0.4 
momentum (A ) 
Parabolic fit to the Cu(1 1 0) surface state 

FIG. 7. Detailed experimental band structure of the Cu(lll) surface 
state. Overlayed, dashed yellow line, is the parabolic fit used to de- 
termine the values in Table [J 

is a uniform doping of the graphene. 

D. Cu(lll) surface state 

An important feature of the graphene on Cu(lll) Fermi 
surface (Fig. [4] (a)) is the well known Cu surface state at 
T indicated by the small circle (dashed yellow line). Previ- 
ous measurements 30 have shown that 1 monolayer of Ag on 
Cu(l 1 1) reduces the bandwidth of this state by a factor of 0.55 
and renormalizes the effective mass by a factor of 0.88. The 
spectral function of the surface state is shown in Fig. [7] with a 
parabolic fit to the state shown in yellow (small dashed lines) 
in Fig. [4] (a) and Fig. [7] Parameters from the parabolic fit are 
used to determine the bandwidth and the Fermi vector pre- 
sented in Table |l| along with the values found for pure Cu and 
1 ML Ag on Cu by Bendounan et al! 3 -^. The bandwidth shift 
is not accompanied by a Fermi velocity renormalization and 
is attributed to a simple charge transfer between the graphene 
andtheCu(lll). 

As both the surface state and the graphene ir band have an 
approximately circular Fermi surface the amount of charge 
transfer can be related to the change in radius of the Fermi 
surface in both cases. Such analysis shows that the amount of 
charge transfer to the surface state is appreciably less than to 
the graphene tt bands. It is therefore expected that the major- 
ity of the charge transferred to the graphene layer comes from 
the Bulk Cu and not the surface layer. 

The presence of a surface superstructure has been shown 
to alter the effective mas d 30 * 31 * as well as open a gap in the 
surface state at the surface brillouin zone boundary ( ~ 0.151 
A -1 for 1 ML Ag on Cu(lllj23). The absence of such a gap 
in the spectral function (Fig. [7]) and the similarity of the ef- 
fective mass from clean Cu(lll) , ~ 0.41, and graphene on 
Cu(lll), ~ 0.39, indicates that no surface superstructure ex- 
ists in this system. This is in agreement with the LEEM data, 
which indicated an oxygen superstructure on the graphene on 
Cu(100) sample, but not on the graphene on Cu(lll) sam- 
ple. It is therefore likely the presence of ordered, intercalated 



TABLE I. Cu(lll) surface state data. The effective mass is un- 
changed by the graphene overlayer, in contrast to 1 ML Ag on 
Cu(lllJ^. 



Substrate 


Bandwidth 
(eV) 


Fermi vector Effective 
(A -1 ) mass 

(m*/irie) 


Cu(lllP 


0.435 


0.215 


0.41 


1ML Ag on Cu(lllpl 


0.241 


0.151 


0.36 


lMLGon Cu(lll) 


0.218 


0.150 


0.39 



oxygen that leads to the second doping level in the graphene 
on Cu(100) samples, which is not observed for graphene on 
Cu(lll). 



IV. CONCLUSION 

We investigated graphene grown on single crystal Cu(lll) 
and Cu(100) surfaces. The graphene layers show a high de- 
gree of rotational disorder resulting in graphene K points 
forming a nearly unbroken arc in the diffraction space mea- 
surements of both ARPES and LEED measurements, this ro- 
tational disorder was also observed in STM measurement^. 
The graphene lattice does preferentially align with the 
Cu(lll) lattice, though. The doping level, measured by the 
offset of the Dirac crossing from the Fermi level (~ -0.3 eV), 
was found to be similar for both substrates and to that found 
for graphene on 1ML Cu/Ni(lll) (~ -0.31 eV). On exposure 
to air oxygen was found to intercalate under the graphene 
on the Cu(100) surface forming a superstructure, and con- 
sequently a second doping level of the overlaying graphene. 
The doping level of this second state is higher (~ -0.6 eV) 
owing to more charge transfer from the oxygen than the Cu. 
The size of the gap induced by interactions with the substrate 
was found to be larger than for graphene on 1ML Cu/Ni(l 11) 
(250 meV and 180 meV respectively) and even larger for the 
intercalated oxygen regions (~ 350 meV). Interestingly the 
interaction between graphene and the Cu(100) and Cu(lll) 
surfaces appears to be similar due to identical doping levels 
and gaps, with a stronger interaction observed with the inter- 
calated oxygen. Similar to previous studies 14 we find that the 
current DFT theory underestimates the band gap (~ 1 1 me\fl3) 
by a factor of 30. The Cu surface state at T on the Cu(lll) 
surface is also found to be protected under the graphene, with 
only a doping derived binding energy shift being observed. 
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